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ABSTRACT 

Although  there  are  sixteen  elements  of  the  Stokes  Matrix,  they  are  con¬ 
structed  from  basically  four  amplitudes  and  three  phase  differences.  This 
of  course  implies  that  there  exists  nine  independent  relationships  connecting 
the  elements.  These  relationships  are  equalities  for  scattering  by  a  single 
particle  in  a  fixed  orientation  and  in  a  fixed  direction.  They  become 
inequalities  when  Stokes  matrices  from  an  ensemble  of  particles  differing 
in  size,  orientation,  morphology,  or  optical  properties  are  added  incoherently. 
These  relations  will  prove  to  be  very  useful  for  providing  consistency  checks 
on  experimental  measurements  of  all  sixteen  elements. 


I.  Introduction 

In  the  classic  work  of  Van  de  Hulst\  he  points  out  the  fact  that  there 

exists  nine  independent  relationships  (equalities)  between  the  sixteen 

elements  of  the  4x4  Stokes  matrix.  These  relations  however  were  not  presented. 

He  also  stated  that  when  Stokes  matrices  were  added  incoherently  that  the 

equalities  became  inequalities.  This  situation  arises  when  one  is  measuring 

the  combined  Stokes  matrix  from  an  ensemble  of  particles  differing  in  size, 

orientation,  morphology  or  optical  properties. 

Experimentalists  are  now  routinely  measuring  aV[  sixteen  elements  of  the 

2-4 

Stokes  matrix  for  both  single  particles  and  collections  of  particles. 

It  is  therefore  important  that  these  relations  be  employed  routinely  to 

check  for  consistency.  To  our  knowledge  the  only  authors  who  have  presented 

5 

the  equalities  explicitly  were  Abhyankar  and  Fymat.  We  will  show  that 
the  first  six  of  their  relations  (10c  -  15c),  which  are  quadratic  in  the 
elements  are  correct,  however  their  remaining  three  equations  (19c  -  21c)  are 
quartic  relations  which  we  will  show  are  actually  the  product  of  two  quad¬ 
ratic  equalities,  one  of  which  is  not  independent  and  is  therefore 
redundant.  We  will  also  present  proofs  for  the  inequalities  along  with 
examples. 

II.  DERIVATION  OF  EQUALITIES 

We  will  use  the  notation  of  Van  de  Hulst1  in  the  following  derivations. 
Let  the  complex  scattered  amplitudes  be  written  as 

2.3,4.  0) 


and  let 


e  - 

81  " 

$2  • 

(2a) 

6  = 

*3  - 

(2b) 

Y  * 

61  - 

64* 

(2c) 

0  = 

84  - 

&2  * 

(2d) 

X  = 

*1  - 

33. 

(2e) 

n  = 

*4  * 

03- 

(2f ) 

With  this  notation  the  elements  of  the  Stokes  matrix  (f^..)  can  be  written 
as  follows: 


fll  =  (“l*  +  a2  *  a3  +  a42^2  * 

(3a) 

f12  =  ('al  +  a2  *  °3  +  “4  )/2  * 

(3b) 

f  1 3  *  «2d3  COS6  +  0^4  COSy  , 

(3c) 

fj4  =  -  a2a3  s*n^“  a|a4  Siny * 

(3d) 

£  t  2.2.2  2\  jm 

f21  =  ('al  +  a2  +  a3  '  °4  )/2  * 

(3e) 

„  ,2.2  2  2\/9 
f22  *  <al  +  °2  -  a3  •  °4  )/2  ‘ 

(3f ) 

^23  =  a20t3  cos8  ‘  ala4  C0SY  » 

(3g) 

^24  =  ~a2a3  Sin5  +  ala4  S^ny  * 

(3h) 

f 31  =  «2a4  COSa  +  ala3  cos*  * 

(31) 

fj2  =  a2a4  COSa  "  ala3  cos*  * 

(3j ) 

fj3  =  a'|ot2  C0Se  +  a3a4  C0Sn  * 

(3k) 

fj4  =-a'ja2  s^ne  +  a3°4  s^nT1 » 

(31) 

^41  *  a2a4  s^na  +  0la3  * 

(3m) 

f 42  -  a2a4  sino  -  <*i«3  sinx  > 

(3n) 

f43  s  al°2  s^ne  +  a3°4  s*nn  * 

(3o) 

(3p) 
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f^4  =  a-|02  COSe  -  0^04  cosn  . 


With  the  elements  written  in  the  above  form,  it  is  quite  straight  forward 
to  obtain  the  following  nine  Independent  equalities: 

^11  +  ^22^  "  ^12  +  ^21  ^  3  ^33  +  ^44^  +  ^43  "  ^34^  =  4otl  “2*  (4a) 

(fH  -  f22)2  "  (f21  “  fi2^  3  ^f33  ‘  f44^  +  ^f43  +  f34^  =  4a32a42 

(fll  +  f21)2  '  (f12  +  f22)2  =  (f13  ♦  f23)2  +  (f14  +  f24)2  =  4a22a32  (4c) 

(fll  f21^  "  ^f12  “  f22^  =  ^f13  "  f23^  +  ^f14  '  f24^  =  4al2a42*  ^4d^ 

(^11  +  ^12^  ”  ^21  +  ^22^  =  ^31  +  ^32^  *  ^41  +  ^42^  *  4a2  a4  ^4e^ 

(fll  "  f12^  ’  ^f21  "  f22^  *  ^f31  "  f32^2  +  ^f41  "  f42^  =  4al2a32  ^4f) 

f13  f14  ‘  f23  f24  3  f33  f34  +  f43  f44  “  "2ala2a3a4  sin  (6}-62+63‘64^  »  ^ 

f14  f23  "  f13  f24  3  f42  f31  ’  f41  f32  3  ■2“l0l2“30l4  sin  ( 6-, +B2-63“ e4 ^  .  (4h) 

f31  f41  ’  f32  f42  3  f33  f43  +  f34  f44  3  2“la2a3a4  sin  (ere2‘e3+B4^  •  (4i) 

In  place  of  eqns.  (4g)  -  (4 i )  we  can  use  the  following  three  independent 

equations,  namely, 

f33?  “  f342  +  f432  ’  f442  “  f13?  "  f142  '  f232  +  f24?  3  4“lot2a3“4  cos  ^l'W^ 

(5a) 

f33?  "  f43?  +  f342  ’  f442  3  f312  "  f41?  ”  f322  *  f42?  =  4ala2a30t4  cos  ^Br82”e3+e4^ 

(5b) 

f312  '  f322  +  f412  ’  f422  3  f142  '  f24?  +  f132  "  f232  =  4ola2ot3a4  cos  ( 6] +B2- 63- 64 ) 

(5c) 

It  should  be  noted  that  eqns.  (4a)  -  (4f)  were  also  obtained  by  Abhyankar  and 
Fymat®;  however,  eqns.  (4g)  -  (4i )  and  (5a)  -  (5c)  were  given  as  products. 

For  example  their  eqn.  (19c)  was  given  as  the  product  of  our  eqns.  (4g)  and 
(5a)  yielding  a  quartic  relation,  which  we  now  see  is  unnecessary.  In  fact. 

It  Is  not  too  difficult  to  show  that  the  square  of  the  lhs  of  eqn.  (5a)  is 
equal  to  the  product  of  the  middle  termsof  eqns.  (4c)  and  (4d)  minus  four 


times  the  square  of  the  middle  term  of  eqn.  (4g). 

An  interesting  equality  can  be  obtained  by  summing  eqns.  (4a)  -  (4f)  namely 

i>:j  f?j  =  4fl/‘  iJ  =  1,2’3‘4-  (6) 

It  is  interesting  to  note  that  for  a  spherical  particle  where 
f31  =  f32  =  f13  =  f23  =  f14  =  f41  =  f42  =  f24  =  0  . 


and 


F11  =  f22*  f12  =  f21’  f33  s  f44’  f34  =  ~f 
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(7a) 

(7b) 


that  only  one  non-trivial  equality  exists  which  is  eqn.  (4a)  which  yields 
'll*  .  f„*  *  f33Z  *  fJ42  (■> 

III.  DERIVATION  OF  INEQUALITIES 

We  now  would  like  to  consider  what  happens  when  we  add  incoherently  the 
Stokes  matrices  from  an  ensemble  of  particles  differing  in  size,  orientation, 
morphology  or  optical  properties.  We  will  use  a  superscript  to  denote 
scattering  amplitudes  and  phases  for  each  member  of  the  ensemble.  With  this 
notation  the  1 hs  (left  hand  side)  of  eqn.  (4a)  becomes 

( f  1  y  +  f22^  ~  ^f12  +  f21^  =  V,'j  “i1  “i1  <l2J  »  (9a) 


whereas  the  rhs  becomes 


(f33  +  f44*2  +  *f43  ‘  Vj'M  C0S  ‘  ' 


(%) 


It  should  be  understood  that  fkj  as  used  here  denotes elements  of  the  total 

Stokes  matrix  of  the  ensemble  and  hence  the  use  of  the  superscript  T  and  also 

that  the  sumnation  indices  i  and  j  are  over  all  members  comprising  the 

ensemble.  Now  the  rhs  of  eqn.  (9a)  can  be  written  as 

-  „  1  i  j  j  .  , 

4  E  a,  a,  a,  a»  -  i 

i,j  1  1  2  2  1 

=  21r.j(o11a2J  -  a^ag1)  +  ^“l^ 

>  4  I  >  4  I  al^a2^al^n2^ 

1 , j  1»j 


p/  i  i  j  j  .  j  j  1  1  \ 

Ela^  <*2  a2  +  °1  al  a2  ”2 


(9c) 


which  is  the  rhs  of  eqn.  (9b);  and,  therefore 

<f,I  *  f,I)2  -  (f,2  *  f2I>2  i  <f33  +  fJ>2  +  (f43T-  f34>2. 


It  is  now  clear  that  the  inequalities  corresponding  to  eqns.  (4a)-(4f) 


become 


*f43  ‘ 

(10a) 

<f43T+ 

f34^. 

(10b) 

<f14T  + 

hP. 

(10c) 

hP, 

(10d) 

<f41T  + 

f421>2. 

(lOe) 

<f4,T- 

f42^. 

OOf) 

Now  the  inequality  for  eqn.  (6)  follows  immediately  and  is 

A  (fij)2±4(fn>2  =  1 ,2,3,4.  (11) 

1  ,j 

One-way  inequalities  for  eqns.  (4g)  -  (4i)  and  (5a)  -  (5c)  cannot  be 
established.  For  example  the  lhs  of  eqn.  (4g)  yields 

2  1  VV«3V  sin  +63J  '  O  * 

*  J 

whereas  the  middle  term  of  eqn.  (4g)  yields 

2  L  a1ia2ia3j«4j  Sin  U-,3  -  ^  -  34’  )^ 

and  we  can  see  that  the  inequality  can  go  in  either  direction  depending  on  the 
amplitudes  and  relative  phase  difference  at  each  scattering  angle  of  particles 
In  the  ensemble. 

It  should  also  be  noted  that  the  inequalities  derived  above  also  hold  if 
one  makes  measurements  on  a  single  particle  but  averages  over  a  finite  solid  angle. 


IV.  Test  Of  Inequalities 

To  validate  the  inequalities  we  have  presented,  we  performed  the 
following  two  calculations.  We  first  computed  the  elements  of  the  Stokes 
matrix  for  a  single  sphere  of  size  parameter  7.9109  with  a  real  refractive 
index  of  1.6146.  The  results  were  then  averaged  over  a  3°  angular  range. 

We  then  tested  eqn.  (10)  by  using  the  normalized  Stokes  matrix  i.e.,  we 
define  ■  -F^ -j  and  then  eqn.  (10)  becomes 

<  4  ;  l.j  -  I.2.3.4.  (14) 

The  results  of  this  calculation  are  shown  in  Fig.  1  along  with  the  only 

non-zero  normalized  elements  (with  the  obvious  exception  of  F,j  *  1  for  all 

scattering  angles)  F12  =  F21 ,  F33  *  F44,  and  F43  =  -F34>  It  is  immediately 

clear  that  there  is  a  great  deal  of  structure  in  the  deviation  of 
2 

.1.  F. .  from  four.  In  fact  a  deviation  as  large  as  18%  occurs  at  a 

1  *  J  I  J 

scattering  angle  of  '  157°. 

The  next  case  we  considered  was  to  average  over  a  polydispersion  of 
spheres  and  test  the  inequality  with  no  angular  averaging.  To  do  this  we 
assumed  a  Gaussian  distribution  of  spherical  particles  with  a  mean  size 
parameter  of  7.9109,  the  same  as  for  the  single  sphere  in  Fig.  1,  and  a 
standard  derivation  of  0.117  in  size  parameter  and  the  same  refractive  index 
as  used  in  the  single  sphere  case.  These  results  are  shown  in  Fig.  2. 

What  is  particularly  noteworthy  is  the  fact  that  even  though  the  normalized 

2 

elements  are  quite  similar  between  the  two  cases,  the  i  F^  •  is  quite 

i  J 

dissimilar  in  structure. 


Due  to  the  quadratic  nature  of  the  inequalities,  we  conjecture 


that  the  use  of  them  may  provide  a  more  discriminating  test  of  small 
differences  between  Stokes  matrix  elements  of  two  slightly  dissimilar 
objects.  Experimental  verification  of  this  idea  is  presently  being 
pursued  and  the  results  will  appear  in  a  future  publication. 

V.  Conclusion 

We  have  presented  both  the  equalities  and  inequalities  which  occur 
between  the  sixteen  elements  of  the  Stokes  matrix.  The  use  of  them  should 
become  a  routine  part  of  any  experiment  designed  to  measure  them.  The 
reason  being  that  they  can  be  used  as  a  self  consistency  test  on  the 
experimental  measurements.  For  example  if  one  finds  that  any  one  of 
the  inequalities  is  violated,  then  one  can  safely  assume  that  there  are 
noise  and/or  systematic  problems  in  the  measurements.  This  appears  far 
superior  than  simply  trying  to  fit  some  theoretical  curve  to  the  data  in 
order  to  determine  the  noise  level.  The  inequalities  also  hold  for  measure¬ 
ments  made  on  a  single  member  of  an  ensemble  averaged  over  a  finite  solid 
angle  which  is  always  realized  experimentally.  We  would  also  like  to 
make  a  conjecture  that  the  inequality  tests  may  provide  a  very  useful  way 
for  discriminating  between  two  slightly  dissimilar  objects. 


This  research  was  supported  by  the  Office  of  Naval  Research  through 
contract  N00014-80-C-0113. 
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FIGURE  CAPTIONS 


Fig.  1 


Fig.  2 


2 

Plot  of  i  F. .  ;i,j  =  1,2, 3, 4,  for  a  single  sphere  as  a  function 
i,j  J 

of  scattering  angle  for  a  3°  angular  average.  The  computations 
were  carried  out  using  Mie  theory  for  a  size  parameter  of  7.9109. 
and  a  real  refractive  index  of  1.6146.  Also  shown  are  the 
normalized  elements  F12,  F33,  and  F43. 


Same  as  Fig.  1  except  a  Gaussian  size  distribution  of  spheres  was 
used  with  a  mean  size  parameter  of  7.9109and  a  standard  deviation 
of  the  size  parameter  of  0.117.  No  angular  averaging  was  done 
for  this  calculation. 
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